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ABSTRACT 

We present the serendipitous discovery of 16 ms pulsed X-ray emission from the 
Crab-like supernova remnant N157B in the Large Magellanic Cloud. This is the fastest 
spinning pulsar associated with a supernova remnant (SNR). Observations with the 
Rossi X-ray Timing Explorer (RXTE), centered on the field containing SN1987A, reveal 
an X-ray pulsar with a narrow pulse profile. Archival ASCA X-ray data confirm this 
detection and locate the pulsar within 1' of the supernova remnant N157B, 14' from 
SN1987A. The pulsar manifests evidence for glitch(es) between the RXTE and ASCA 
observations which span 3.5 years; the mean linear spin-down rate is P = 5.126 x 
1(C 14 s s _1 . The background subtracted pulsed emission is similar to other Crab-like 
pulsars with a power law of photon index of ~ 1.6. The characteristic spin-down age 
(~ 5000 years) is consistent with the previous age estimate of the SNR. The inferred 
.B-field for a rotationally powered pulsar is ~ 1 x 10 12 Gauss. Our result confirms the 
Crab-like nature of N157B; the pulsar is likely associated with a compact X-ray source 
revealed by ROSAT HRI observations. 

Subject headings: pulsars: general — pulsars: individual (PSR J0537— 6910) — X-rays: 
general — supernova remnant 
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1. Introduction 

Crab-like supernova remnants (SNRs) play a 
critical role in our understanding of young pul- 
sars and their pulsar wind nebulae. These rare 
SNRs which contain central pulsars are distin- 
guished by their centrally-filled morphologies and 
non-thermal X-ray spectra. The X-ray emission 
of such SNRs comes predominantly from a syn- 
chrotron nebula powered by an embedded young 
pulsar (e.g., Seward 1989). Only three confirmed 
members of this class were previously known: the 
Crab Nebula with its famous 33 ms pulsar PSR 
B0531+21, the SNR B0540-693 in the Large 
Magellanic Cloud (LMC) with its 50 ms pul- 
sar PSR B0540— 69, and the nebula around the 
150 ms pulsar PSR B1509— 58. Several candi- 
date Crab-like SNR have been reported. These 
have similar morphologies and spectral charac- 
teristics but so far lack a detected pulsar. Young 
Crab-like SNRs are expected to evolve into com- 
posite SNRs, in which the thermal emission from 
shock-heated gas will rival or exceed the dimin- 
ishing X-ray radiation from the pulsar-powered 
synchrotron nebulae. This evolutionary sequence 
is not yet firmly established, though, mainly be- 
cause of the limited number of samples. 

We report in this Letter the discovery of an 
ultra fast pulsar in N157B (also known as NGC 
2060, SNR 0538-69.1, & 30 Dor B; Henize 1956), 
confirming the Crab-like nature of this remnant 
(Wang & Gotthelf 1998 and ref. therein). The 
pulsar was first detected serendipitously with 
RXTE data while searching for a pulsed signal 
from nearby SN1987A (Marshall et al. 1998). We 
used archival ASCA data of the region to confirm 
the detection of the pulsed emission, to measure 
its spin-down rate, and to locate the emission to 
within the boundary of N157B. Based on the 
ASCA position, we will refer to the pulsar as 
PSR J0537— 6910. In the following, we present 
our discovery of PSR J0537— 6910 and discuss its 
properties in the context of other young pulsars. 



2. Observations 

The LMC field containing SN1987A has been 
frequently sampled by many X-ray observatories 
including RXTE and ASCA. There were two ob- 
serving runs in 1996 with RXTE to search for 
pulsations from SN1987A (see Observation Log, 
Table 1). The first was an uninterrupted obser- 
vation of 21 ks on 12 Oct. 1996. The other run, 
which began on 20 Dec. 1996, obtained a total 
of about 78 ks of good observing time during an 
elapsed time of 184 ks. 

The RXTE Observatory (Bradt et al. 1993) 
consists of five co-aligned collimated detectors 
known collectively as the Proportional Counter 
Array (PC A), a set of crystal scintillation de- 
tectors known as the High Energy X-ray Timing 
Experiment (HEXTE), and the All-Sky Moni- 
tor. Here we report exclusively on data obtained 
with the PCA because of its large effective area 
(~ 6500 cm 2 at 10 keV). The PCA (Jahoda et al. 
1996) has collimators that produce a roughly cir- 
cular aperture with a ~ 1° FWHM response. 
There are several known X-ray sources within 
the field-of-view (FOV) with flux above 2 keV, 
the brightest of which are the black-hole candi- 
date LMC X-l, the 50 ms pulsar PSR B0540-69, 
and N157B, the object of this study. The PCA 
observations used the Good Xenon data mode, 
which time-tags each photon with 0.9/Us timing 
resolution. For the current analysis, the abso- 
lute timing uncertainty is ~ 100/xs (Rots et al. 
1998). Moderate spectral information is available 
in the 2 — 60 keV energy band with a resolution 
of ~ 16% at 6 keV. 

The ASCA Observatory (Tanaka et al. 1994) 
acquired three long ~ 20 ks observations of the 
region around SN1987A with sufficient timing 
resolution to detect ms pulsations (see Table 1). 
These were all obtained in 1993 during the early 
Performance and Verification phase of the ASCA 
mission. Two of the observations containing 
N157B were optimal, placing the SNR close to 
the optical axis, while the third data set observed 
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the SNR with reduced efficiency, at 18' off-axis 
where the source flux is strongly vignetted. 

ASCA data were acquired with all four co- 
aligned focal plane instruments. We used data 
exclusively from the two Gas scintillation Imag- 
ing Spectrometers (GIS), and only data for which 
time-tagged photons were available in the high- 
est time resolution mode (61 [xs or 0.488 ms de- 
pending on telemetry rate). ASCA GIS timing 
measurements are found to have an absolute ac- 
curacy of 200 [is in this mode (see Saito et al. 
1997). Data from the Solid State Imaging Spec- 
trometers (SIS) do not have adequate time reso- 
lution for the current study. The GIS offers ~ 8% 
spectral resolution at 6 keV over its ~ 1 — 12 keV 
energy band-pass. Each GIS sits at the focus of 
a conical foil mirror and the combination results 
in a spatial resolution of 1 — 3 arcmin (depending 
on energy) over the GIS's ~ 50' arcmin diameter 
active FOV. Point sources can be located with 
~ 0.6' accuracy (see Gotthelf 1994). 

Both the RXTE and ASCA data were edited 
to exclude times of high background contami- 
nation using the standard processing screening 
criteria ( "REV2" for ASCA). This rejects time 
intervals of South Atlantic Anomaly passages, 
Earth occultations, bright Earth limb in the field- 
of-view (ASCA only), and other periods of high 
particle activity. For each observation, event 
data from all detectors were co-added and the 
arrival times of each photon corrected to the 
solar system barycenter using the JPL DE200 
ephemeris. 

3. Timing Results 

For this Letter we have concentrated on RXTE 
data collected in the 2 — 10 keV energy band in 
layer 1 of the PCA, which provides the best sen- 
sitivity for weak sources with a Crab-like spec- 
trum, and the ~ 10 — 25 keV hard energy band 
using all three PCA layers. We also examined 
data above 30 keV, which are expected to be 
dominated by detector background, to verify that 



there are no similar periodicities produced by the 
instrument itself. 

We searched for significant signals using 2 27 
point Fast Fourier transforms on data binned in 
0.16 ms steps. The 12 Oct 1996 observation pro- 
duced a 7.5a signal at ~ 62.055 Hz in the hard 
energy band. Similar power was found for the 
2 nd and 3 rd harmonics. The probability of false 
detection is < 1 in 10 10 . This result is repro- 
duced in the 2 — 10 keV band data, with a IOct 
detection in the first three harmonics of the 62 
Hz signal. No other significant power was found 
with the exception of frequencies corresponding 
to the 50 ms pulsar PSR B0540-69. The new 
frequencies are not simple harmonics or aliases 
of the 50 ms pulsar signal. As we show below, 
any association with the 50 ms pulsations can be 
ruled out conclusively. We folded the data into 20 
bins on periods near 16 ms and computed \ 2 f° r 
a model with a count rate independent of pulse 
phase. The largest value for \ 2 of 157 (7.8a) 
was found for a period of P = 16.114717 ms. 
The pulsar was also detected during the second 
RXTE observing run, allowing us to derive an 
initial period derivative. 

We examined the ASCA data for the RXTE 
pulsations. N157B was considered the most likely 
source of the pulsations because it is the bright- 
est source in the RXTE FOV after LMC X-l and 
B0540— 693. It is a barely resolved, moderately 
bright ASCA source in the 30 Doradus star for- 
mation region. N157B is clearly imaged in the 
ASCA GIS as an isolated source, ~ 10' south- 
west of the X-ray bright R136 region, and ~ 15' 
northwest of the 50 ms pulsar PSR B0540— 69, 
which produced stray light in the ASCA FOV 
containing N157B. 

We extracted photons from both GIS instru- 
ments using a 4' radius aperture centered on the 
peak flux from N157B and constructed a peri- 
odogram as with the RXTE data. For each trial 
period, we folded the data into 10 bins to opti- 
mize the signal in the photon noise dominated 
data and computed the x 2 of the resultant pro- 
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file. From each ASCA observation we conducted 
a high resolution blind search of the data for sig- 
nificant power about a range of ±0.5 ms around 
the expected period. To search for short duty cy- 
cle pulsations in period-space we oversampled the 
period resolution by a factor of 20, using incre- 
ments of 0.05 x P 2 /T, were T is the observation 
duration, and P is the test period. 

A highly significant signal (> 12a) was found 
during both on-axis ASCA observations. The low 
signal-to-noise September 1993 observation pro- 
duced a 5a result, whose reduced significance is 
consistent with the target being placed 18' off- 
axis. This is the only significant peak found from 
all data sets, over the range searched. We used 
the change in the period between the two on- 
axis observations to bootstrap the period deriva- 
tive, which resulted in a search of greater sen- 
sitivity. We also find that the period detection 
is greatly improved by excluding photons below 
< 2 keV. This suggests that the pulsar signal 
is significantly contaminated by un-pulsed, soft 
thermal flux from the SNR nebula, which domi- 
nates below 2 keV. The non-detection of the pul- 
sar by Wang & Gotthelf (1998) is explained by 
the narrowness of the peak, the duration of the 
observations (typically 90 ks), and the lack of a 
period derivative a priori. 

Since it is not possible to maintain phase be- 
tween any of the four observation intervals, esti- 
mates of the period derivative are based on the 
change in period between observations. Figure 1 
shows the period evolution and the residuals from 
the best-fit linear model of P = 0.016113089 s ± 
1.7 ns at MJD 50000 (at Earth) and P = 5.126± 
0.002 x 10~ 14 s s -1 . Our estimates of the pe- 
riod uncertainty for each data point in Figure 1 
are derived by Monte-Carlo simulation and Max- 
imum Likelihood analysis. Clearly, there are sig- 
nificant deviations from a constant period deriva- 
tive; these are not removed by including a pe- 
riod second derivative to the fit. This inability to 
find a single solution for both sets of observations 
suggest that one or more glitches in the period 



has occurred (see discussion in §4). The mea- 
sured period derivatives during the two epochs 
are Pasca = 5.154 ± 0.002 x 10~ 14 s s" 1 and 
Phxte = 5.174 ± 0.002 x 10~ 14 s s _1 . 

Using the above ephemeris, we generated pulse 
profiles for the 1993 ASCA and the 1996 RXTE 
combined data sets. The resulting profiles are 
displayed in Figure 2. They are well character- 
ized by a single, narrow peaked pulse with an ap- 
proximately Gaussian shape of FWHM of ~ 1.7 
ms (10% duty cycle) ms. The profiles appear un- 
changed between the 1993 and 1996 observations. 
The pulsed emission comprises 9.6% of the total 
ASCA counts in the 4' radius aperture above 2 
keV. 

To demonstrate that these pulsations are uniquely 
associated with N157B, an image of the pulsed 
emission was generated by subtracting the off- 
pulsed data from the on-pulsed image. Only 
a single source of significant emission remained, 
the pulsed emission, an unresolved ASCA point- 
like source located at the coordinates of the SNR 
N157B. We thus reproduce and verify the RXTE 
measurement, and unambiguously locate the ori- 
gin of the pulsed emission to within an arc minute 
of the region encompassing N157B. 

In summary, we consider this a conclusive de- 
tection of a new pulsar in N157B, PSR J0537-6910, 
whose period of 16 ms makes it the fastest known 
pulsar associated with a SNR. 

4. Discussion 

The steady increase in the period suggests 
that the pulsed X-ray emission is powered by the 
spin-down energy loss of the pulsar. We esti- 
mate the spin-down power as E = (2ir) 2 IP/P 3 = 
4.8 x 10 38 /45 er g s s" 1 ) where Z45 is the moment of 
inertia of the neutron star in units of 10 45 gcm 2 . 
For a pulsar losing rotational energy via mag- 
netic dipole radiation, E ~ (B 2 R 6 Q 4 )/(6c 3 ), and 
the surface magnetic field strength at the pole is 
B p ~ 1.0 x 10 12 G. Here, R ~ 10 km is the neu- 
tron star radius, Q = 2ir/P is the angular veloc- 
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ity of the rotation, and c is the speed of light. 

The pulsed spectral component can be isolated 
using phase-resolved spectroscopy. The source 
photons were folded at the best measured period 
and phase dependent spectra constructed. Us- 
ing the off-pulse spectrum as background, we fit 
the on-pulsed spectrum with a simple absorbed 
power law model to the ASCA &i RXTE data in 
the 2 — 20 keV band. The best fit photon index 
(xl = 0.96 for 30 DoF) is 1.6 (1.3-2.0), which is 
consistent with the typical value ~ 1.7 for other 
pulsars (Lyne & Graham-Smith 1990). The large 
uncertainties are due to subtracting two spec- 
tra with similar number of counts. The unab- 
sorbed 2-10 keV pulsed flux of PSR J0537-6910 
is ~ 6.7 ± 0.6 x 10 -13 ergs s _1 cm -2 (see Ta- 
ble 1). The corresponding pulsed luminosity is 
~ 1.7 x 10 35 ergs s _1 (into Air) assuming a dis- 
tance of 47 kpc (Gould 1995). 

The derived pulsed luminosity suggests that 
~ 4 x 10~ 4 of the pulsar's spin down energy is 
emitted as pulsed X-rays in the 2 — 10 keV band. 
This luminosity is consistent with the empiri- 
cal relation, L x « 7 x W 26 (B 12 /P 2 ) 2 - 7 ergs s -1 , 
which fits well to the magneto-rotation driven 
X-ray emission from previously known pulsars 
(Ogelman 1995), to within the uncertainty in the 
power index of 2.7 ± 0.5. Here, B\2 is the mag- 
netic field in units of 10 12 G and P 2 is the spin 
period in units of seconds. 

Compared to other young pulsars known to be 
associated with Crab-like SNRs, PSR J0537-6910 
is quite unusual. Its X-ray pulse profile appears 
to be the narrowest among the young pulsars. 
The Crab pulsar is double peaked and the profiles 
of both PSR B0540-69 and PSR B1509-58 are 
very broad (duty cycle > 0.3). The pulse width 
of PSR J0537-6910, 1.7 ms, indicates that the 
size of the emission region is likely smaller than 
about half of the light cylinder radius. A detailed 
study of the pulse profile may help to place con- 
straints on pulsar emission models. The charac- 
teristic age of the pulsar is r = Pf (2P) = 5 x 10 3 
yrs. This age is much greater than any of the 



other Crab-like pulsar, yet similar to the age es- 
timate for N157B of ~ 5 x 10 3 yrs, based on 
X-ray measurements of the size and temperature 
of the remnant (Wang & Gotthelf 1998). An up- 
per limit on the age of the remnant, based on the 
kinematics of Ha-emitting gas in the region (Chu 
et al. 1992), suggests that the remnant may be 
as old as ~ 2 x 10 4 yrs. 

An accurate determination of the remnant age 
is critical to estimate the initial rotation rate of 
the pulsar. PSR J0537— 6910 rotates twice as 
fast as the Crab pulsar, and was likely spinning 
much more rapidly when it was born, depending 
on its assumed braking index and age. Following 
Kaspi et al. (1997), we plot in Figure 3 the possi- 
ble initial spin period vs. age for a range of brak- 
ing index. If the older age estimate is correct and 
if glitches are not important in the overall period 
evolution of the pulsar, the braking index would 
then have to be unusually small (< 1.2) for an 
initial period > 1 ms. If the younger value for the 
age is more likely, the initial spin is then a few ms, 
assuming n ~ 3, as in the case of the magnetic 
dipole model and for other young pulsars. In 
either case, these values provide important con- 
straints on neutron star birth spin models. 

The phenomena of pulsar glitches is well es- 
tablished and provides important constrains on 
the moment of inertia of neutron stars (see Shapiro 
& Teukolsky 1983). These sudden changes in 
the pulsar spin period and period derivative are 
attributed to "starquakes" , stress relief between 
the neutron star crust and the superfluid core. 
Monitoring of the Crab and Vela pulsars, as well 
as other young pulsars, show that large glitches 
occur on ~ a few years timescales. There is com- 
pelling evidence that PSR J0537-6910 has un- 
dergone one or more glitches between the ASCA 
and RXTE epochs (see figure 1). The integrated 
magnitude of the glitches, Af2/f2 ~ 10 -6 , is com- 
parable to that observed from the Vela pulsar, 
but substantially larger than any of those ob- 
served from the Crab pulsar. More careful analy- 
sis of the data is still needed, however, to quantify 
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the presence of any glitches and to the distinguish 
them from timing noise. 

The origin of the un-pulsed X-rays detected 
from N157B in the energy band > 2 keV is likely 
a synchrotron nebula powered by the pulsar via a 
relativistic wind (Gallant k Arons 1994 and ref- 
erences therein). Recent work by Wang k Got- 
thelf (1998) has shown that N157B contains a 
bright, elongated, and non-thermal X-ray feature 
whose origin is still uncertain. Attached to this 
feature is a compact source with a spatial extent 
of ~ 7", which they claim most likely represents 
the pulsar and its interaction with surrounding 
medium. The X-ray spectrum of N157B, which 
is dominated by the un-pulsed emission, can be 
characterized by a power law with an energy 
slope ~ 1.5 (Wang k Gotthelf 1998), signifi- 
cantly steeper than those (~ 1.0) of other Crab- 
like SNRs (e.g., Asaoka k Koyama 1990). The 
remnant has probably evolved to a stage close 
to that of a composite SNR. This steep spec- 
trum and the elongated morphology suggest that 
something unusual may be happening in N157B. 
Most likely, relativistic particles are transported 
from the pulsar to a radio-emitting wind bubble 
in a collimated outflow (Wang k Gotthelf 1998). 
A substantial fraction of high energy particles 
from the pulsar lose their energy radiatively in 
a bow shock, rather than adiabatically through 
diffusion. 

5. Conclusion 

We have found 16 ms pulsed X-ray emission 
from the 30 Doradus region of the LMC. We 
have shown that these pulsations are associated 
uniquely with the X-ray emission from the SNR 
N157B. This remnant is an important new addi- 
tion to the class of Crab-like SNRs. It provides 
a rare laboratory for studying both an unusually 
rapidly-rotating pulsar and its relativistic wind, 
as well as the structure and evolution of a neu- 
tron star. Future observations are important to 
constrain the age of the pulsar, the period second 
derivative, and to measure putative glitches. 
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Table 1 

Observation Log for PSR J0537-6910 



Data Set 


Date a 


Epoch" 


Exposure/Duration 6 


Period 6 


Pulsed Flux c 




(UT) 


(MJD) 


(ks) 


(ms) 


(xl0~ 13 ergs s _1 cm~ 2 ) 


ASCA 


13 Jun 1993 


49151.306618 


30.8/89.4 


16.1093291 


6.4 ±0.8 


ASCA 


20 Aug 1993 


49219.986830 


32.1/94.8 


16.1096349 


7.1 ±0.9 


ASCA 


23 Sep 1993 


49253.982891 


35.6/82.8 


16.1097864 


4.1 ±2.0 


RXTE 


12 Oct 1996 


50368.148905 


21.8/21.8 


16.1147174 


6.8 ± 1.2 


RXTE 


20 Dec 1996 


50437.550201 


78.0/184.5 


16.1150276 


6.7±0.7 



a For the start of the observation. 

b The 90% confidence level uncertainty on all period measurements is ~ 0.1 ns, except for the 20 Dec 1996 
measurement, which is 0.04 ns. 

c The unabsorbed pulsed RXTE and ASCA flux using an absorbed power law model of photon index 1.6. 
The ASCA GIS flux are derived using photons extracted from 4' radius aperture centered on the ASCA N157B 
position. The off-pulse spectrum was used as a background. See text for details. 
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Fig. 1. — The pulse period evolution of 
PSR J0537-6910 and its residuals from the best 
fit model (see text) assuming a linear spin-down 
trend. 
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Fig. 3. — Predicted initial period 

of PSR J0537— 6910 as a function of the age of 
the pulsar. Here we have assumed a power-law 
deceleration model P oc p 2 ~ n (e.g., Shapiro & 
Teukolsky 1983). The four curves corresponds to 
different values of the index: 1.2 (dot-dashed), 2 
(dotted), 3 (solid), and 4 (dashed). The vertical 
lines indicate a pulsar age of 5 x 10 3 & 20 x 10 3 
yrs old. 



Fig. 2.— The pulse profile of PSR J0537-6910 
in the 2-10 keV band from the RXTE (top) and 
ASCA (bottom) observations. The two profiles 
have been aligned to place the peak emissions at 
the 0.5 phase bin. The relative phases of the two 
measurements are arbitrary. The RXTE profile 
includes 100 ks of data; the ASCA profile is re- 
stricted to the 37 ks of 64 fj,s resolution data. 
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